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Key words: clusters: general - clusters: cosmology 
Abstract. We present a measure of the power spectrum 
on scales from 15 to 800/i _1 Mpc using the ROSAT-ESO 
Flux-Limited X-Ray (REFLEX) galaxy cluster catalogue. 
The REFLEX survey provides a sample of the 452 X-ray 
brightest southern clusters of galaxies with the nominal 
flux limit S = 3.0 x 10~ 12 ergs" 1 cm" 2 for the ROSAT 
energy band (0.1 — 2.4) keV. Several tests are performed 
showing no significant incompletenesses of the REFLEX 
clusters with X-ray luminosities brighter than 10 43 ergs -1 
up to scales of about 800/i -1 Mpc. They also indicate 
that cosmic variance might be more important than pre- 
vious studies suggest. We regard this as a warning not to 
draw general cosmological conclusions from cluster sam- 
ples with a size smaller than REFLEX. Power spectra, 
P{k), of comoving cluster number densities are estimated 
for flux- and volume-limited subsamples. The most impor- 
tant result is the detection of a broad maximum within the 
comoving wavenumber range 0.022 < k < 0.030 /iMpc -1 . 
The data suggest an increase of the power spectral am- 
plitude with X-ray luminosity. Compared to optically se- 
lected cluster samples the REFLEX P{k) is flatter for 
wavenumbers k < 0.05 h Mpc -1 thus shifting the max- 
imum of P(k) to larger scales. The smooth maximum 
is not consistent with the narrow peak detected at k = 
0.05 /iMpc -1 using the Abell/ACO richness > data. In 
the range 0.02 < k < 0.4/iMpc" 1 general agreement is 
found between the slope of the REFLEX P(k) and those 
obtained with optically selected galaxies. A semi-analytic 
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description of the biased nonlinear power spectrum in red- 
shift space gives the best agreement for low-density Cold 
Dark Matter models with or without a cosmological con- 
stant. 



1. Introduction 

The fluctuation power spectrum, P(k), of the comoving 
density contrast, S(r), is a powerful summary statistic 
to explore the second-order clustering properties of cos- 
mic structures. Its direct relation to theoretical quantities 
makes it an ideal tool for the discrimination between dif- 
ferent scenarios of cosmic structure formation and cosmo- 
logical models in general. However, measurements give the 
spatial distribution of 'light' and not the fluctuations of 
the underlying matter field. For galaxies the connection 
between mass and the presence of a stellar system is com- 
plicated because nonlinear gravitational, dissipative, and 
radiative processes could lead to a nonlinear biasing up 
to rather large scales (e.g., Bertschinger et al. 1997 and 
references given therein) . For rich clusters the relation be- 
tween mass and the presence of such systems is expected 
to be governed by comparatively simple biasing schemes 
(e.g., Kaiser 1984, Bardeen et al. 1986, Mo & White 1996), 
mainly driven by gravitation, and only slightly modified 
by dissipative processes. In this sense rich clusters of galax- 
ies are much easier to model and thus 'better' tracers of 
the large-scale distribution of matter. 

Power spectra obtained from optically selected clus- 
ter surveys (Peacock & West 1992; Einasto et al. 1993; 
Jing & Valdarnini 1993; Einasto et al. 1997; Retzlaff et 
al. 1998; Tadros, Efstathiou & Dalton 1998) are found to 
have slopes of about —1.8 for k > 0.05 h Mpc -1 and a 
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turnover or some indications for a turnover at k ~ 0.03 — 
0.05 ZiMpcT 1 . Contrary to this, Miller & Batuski (2000) 
find no indication of a turnover in the distribution of Abcll 
richness > 1 clusters for k > 0.009 h Mpc _1 .[] Measure- 
ments on scales > 500/i~ 1 Mpc or k < 0.013 /iMpc -1 
where the cluster fluctuation signal is expected to be 
smaller than 1 percent are, however, extremely sensitive 
to errors in the sample selection. The resulting artificial 
fluctuations increase the measured power spectral densi- 
ties and thus prevent any detection of a decreasing P(k) 
on these large scales. 

The current situation regarding the detection and the 
location of a turnover in the cluster power spectra ap- 
pears to be very controversal with partially contradict- 
ing results. Physically, the scale of the expected turnover 
is closely linked to the horizon scale at matter-radiation 
equality. This introduces a specific scale into an other- 
wise almost scale-invariant primordial power spectrum 
and thus helps to discriminate between the different sce- 
narios of cosmic structure formation discussed today. The 
narrow peak found for Abell/ACO clusters by Einasto et 
al. (1997) and Retzlaff et al. (1998) suggests a period- 
icity in the cluster distribution on scales of 120/i _1 Mpc 
and, if representative for the whole cluster population, is 
very difficult to reconcile with current structure formation 
models. The undoubted identification of the location and 
shape of this important spectral feature must, however, in- 
clude a clear documentation of the quality of the sample 
from which it was derived. 

Although the quality of optically selected large-area 
cluster samples has been improved during the past years 
by the introduction of, e.g., automatic cluster searches 
(e.g., Dalton et al. 1992, Lumsden et al. 1992, Collins et al. 
1995) a major step towards precise fluctuation measure- 
ments on very large scales is offered by the use of X-ray 
selected cluster samples where also poor systems can be 
reliably identified and characterized within the global net- 
work of filaments or other large-scale structures. This is 
due to several facts. 

First, the relation between X-ray luminosity and total 
cluster mass as observed (see eq.|l^, Reiprich & Bohringer 
1999, Borgani & Guzzo 2000) and as indicated to first or- 
der from the modeling of clusters as a homologous group 
of objects scaling with mass (Kaiser 1986), convincingly 
demonstrates the possibility to select clusters basically by 
their mass, although the la scatter for the determina- 
tion of the gravitational mass from X-ray luminosity is 
still quite large (about 50 percent). This is clearly prefer- 
able compared to a selection of clusters by their optical 
richness, as indicated for example by the results obtained 
within the ENACS (Katgert et al. 1996) where about 10 

1 The Hubble constant Ho is given in units of h = 
i?o/(100kms -1 Mpc" 1 ) and the X-ray source properties (lu- 
minosities, etc.) for h = 0.5, the cosmic density parameter is 
fio = 1, and the normalized cosmological constant Qa = 0. 



percent of the Abell, Corwin & Olowin (1989) clusters 
with z < 0.1 (located in the southern hemisphere) do not 
show any significant concentration along the redshift di- 
rection and must thus be regarded as spurious. 

Second, although the spatial galaxy number density 
profiles are more concentrated towards the cluster cen- 
tres compared to the gas density profiles, it is the much 
more centrally peaked X-ray emissivity profile (~ Pg as ) 
which increases the contrast to the background distribu- 
tion and enhances the angular resolution of an X-ray clus- 
ter survey. This decreases the probability of 'projection 
effects' known to contaminate, e.g., the optically selected 
Abell/ACO cluster sample (Lucey 1983, Sutherland 1988, 
Dekel et al. 1989). 

Third, the large-scale variation of galactic extinction 
modifies the local sensitivity of cluster detection (for the 
optical passband see Nichol & Connolly 1996). In addi- 
tion to galactic obscuration galaxies can be confused with 
faint stars which reduces the contrast of a cluster above 
the background so that the system appears less rich (Post- 
man, Geller & Huchra 1986). The resulting artificial dis- 
tortions must be reduced because they easily dominate 
any measured fluctuation on large scales (e.g., Vogeley 
1998). In the following it will be shown that in X-rays the 
local survey sensitivity can be readily computed using the 
local exposure time of the X-ray satellite and the local 
column density of neutral galactic hydrogen, JVhi. 

First results of a power spectrum analysis using X- 
ray selected subsamples of the 291 clusters of the ROSAT 
Bright Survey (Schwope et al. 2000) are presented in Ret- 
zlaff (1999) and Retzlaff k Hasinger (2000). For the count 
rate-limited subsample indications for a turnover of P(k) 
at k — 0.05/i _1 Mpc are found. For the volume-limited 
subsample the statistical significance of this specific fea- 
ture is very weak or almost absent. 

In this paper we present the results of a power spec- 
trum analysis obtained with a sample of 452 ROSAT ESO 
Flux-Limited (REFLEX) clusters of galaxies. A related 
study of the large-scale distribution of REFLEX clusters 
using the spatial two-point correlation function can be 
found in Collins et al. (2000). Sect.|2|gives a brief overview 
of the selection of the cluster sample. Sect.|| concentrates 
on the discussion of the overall completeness of the RE- 
FLEX sample, drawing special attention to those selection 
effects which might limit the fluctuation measurements on 
large scales. In Sect. || standard methods of power spectral 
analyses are applied to estimate P(k). The systematic and 
random errors are computed using a set of N-body simula- 
tions of an open Cold Dark Matter (OCDM) model which 
is shown to give a good though not optimal representation 
of the REFLEX sample (Sect.g). The results arc shown in 
Sect. U and compared with optically selected cluster and 
galaxy samples. In Sect.|7] a semi-analytic model is de- 
rived and compared with the observed power spectra of 
flux- and of volume- limited subsamples. Sect.|^ summa- 
rizes and discusses the main results. 
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Fig. 1. Spatial distribution of the REFLEX clusters of galaxies. Radial axes are given in units of [kms -1 ]. Plotted are all 
clusters with z < 0.35. Galactic extinction partially obscures the regions 6 h — 10 h and 16 h — 20 h . Note the Shapley concentration 
at R.A. = 13 h and cz w 14000 kms -1 surrounded by a filament of galaxy clusters. 



2. The REFLEX cluster sample 

In the following a brief overview of the sample construc- 
tion is given. A detailed description of the various reduc- 
tion steps, the resulting sample sizes, the methods for the 
X-ray flux, S, and luminosity, Lx, computations, the de- 
termination of temperature- and redshift-dependcnt flux 
corrections, as well as the correlation with optical galaxy 
catalogues, and the computation of the local survey flux 
limits (survey sensitivity) can be found in Bohringer et al. 
(2000a,b). 

The REFLEX clusters are detected in the ROSAT All- 
Sky Survey (Trumper 1993, Voges et al. 1999). They are 
distributed over an area of 4.24 sr (13 924dcg 2 ) in the 
southern hemisphere below +2.5 deg Declination. To re- 
duce incompleteness caused by galactic obscuration and 
crowded stellar fields the sample excludes the area ±20 deg 
around the galactic plane and 0.0987 sr at the Small 
and the Large Magellanic Clouds, basically following the 
boundaries of the corresponding UK Schmidt plates (e.g., 
Heydon-Dumbleton, Collins & MacGillivray 1989). 

The sample is based on an MPE internal source cat- 
alogue extracted with a detection likelihood > 7 from 
the ROSAT All-Sky Survey (RASSII). 54 076 southern 
sources have been re-analysed with the growth curve anal- 
ysis method (Bohringer et al. 2000b) which is especially 
suited to the processing of extended sources. Although 
the data were analysed in all three ROSAT energy bands 
most weight is given to the hard band (0.5-2.0 keV) where 
60 to 100 percent of the cluster emission is detected, the 
soft X-ray background is reduced by a factor of approx- 
imately 4, and the contamination through the majority 
of RASSII sources is lowest, so that the signal-to-noise 
for the detection of clusters is highest. As expected the 
new count rates are systematically higher (up to an order 



of magnitude) compared to the count rates given by the 
standard ROSAT analysis software which is optimized for 
the processing of point-like sources. 

The low source counts of many RASS sources as well 
as the limited spectral resolution of the PSPC do not 
give enough information for a proper identification of 
the sources based only on the X-ray properties so that 
additional reduction steps are necessary. Optical cluster 
counterparts are found using counts of COSMOS galaxies 
(Heydon-Dumbleton et al. 1989) in concentric rings with 
different apertures centered around the X-ray source po- 
sitions. The probability thresholds used for the different 
rings are set low to select also weak excesses of galaxy 
surface number densities above background, introducing 
a formal sample incompleteness of less than 10 percent. 

The cluster candidates are screened using the X-ray, 
optical, and literature data. Obvious multiple detections, 
and candidates with a strong point-like contamination 
(e.g., active galactic nuclei AGN) of the X-ray flux where 
the residual flux from the cluster is estimated to be 
smaller than the nominal REFLEX flux limit, are re- 
moved. Double sources are deblended, and count rates 
measured in the hard band are converted to unabsorbed 
fluxes in the ROSAT band (0.1 - 2.4) keV using standard 
radiation codes for a thermal spectrum with temperature 
k^T = 5.0 keV, rcdshift z = 0, metal abundance 0.3 so- 
lar units, and local A^hi (Dickey & Lockman 1990, Stark 
et al. 1992). The internal errors of the measured fluxes 
range between 10 and 20 percent. The effects of a possi- 
ble systematic underestimation of the total fluxes, mainly 
caused by the incomplete sampling of the outer parts of 
the cluster X-ray emission, are presently investigated (H. 
Bohringer et al., in preparation). For the present investi- 
gation the measured fluxes (not the total fluxes) are used. 
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Fig. 2. Number of REFLEX clusters of galaxies as a function 
of galactic longitude in and latitude bu (steps) compared with 
the number of clusters expected for a random realisation of the 
REFLEX selection function (area and sensitivity), computed 
for the nominal flux limit 3 x 10~ 12 ergs -1 cm~ 2 and the mini- 
mum of 10 source counts (continuous lines) . Narrow count bins 
are chosen to show the effects of large-scale clustering. 



A complete identification of all cluster candidates and 
a measure of their redshifts has been performed in the 
framework of an ESO Key Programme (Bohringer et al. 
1998, Guzzo et al. 1999). During this campaign, 431 X-ray 
targets were observed with an average of about 5 spectra 
per target. 

The iterative computation of the X-ray luminosity uses 
in the first step the redshift and the unabsorbed X-ray 
flux to give a first estimate of Lx- This luminosity and 
the luminosity-temperature relation of Markevitch (1998, 
without correction for cooling flows) is used to improve the 
initial temperature estimate (5 keV) . In the next step the 
count rate-flux conversion factor is recomputed including 
now the effects of z. The cluster restframe luminosity is 
calculated by taking into account the equivalent to the 
cosmic K-correction. The X-ray luminosities are given for 
the (0.1 — 2.4) keV energy band (h — 0.5). For this band 
and for clusters with redshifts z < 0.3 and the temperature 
T = 5keV the K-corrections are less than 12 percent. 
Note that the iterative calculation does not introduce any 
uncertainty in the selection function, since each value of 
Lx has a unique correspondence to the first calculated 



unabsorbed flux and thus to a uniquely determined survey 
volume. 

Adding to the above mentioned selection criteria 
the nominal flux limit of the REFLEX sample, 3 x 
10~ 12 erg s _1 cm~ 2 within the ROSAT energy band (0.1 — 
2.4) keV, we find 452 clusters. Of these 449 have mea- 
sured redshifts, 1 object is clearly a cluster while 2 are 
unconfirmed candidates. 65 percent of the sample are 
Abell/ACO/Supplement clusters. However, note the diffi- 
culty to compare X-ray flux-limited and richness-limited 
cluster samples (see Bohringer et al. 2000a for more de- 
tails). 81 percent of these clusters show a significant X- 
ray extent (determined with the growth curve analysis 
method) . This shows how a selection based solely on X-ray 
extent would have missed, given the quality of the RASS II 
data, a significant percentage of true clusters. Less than 
10 percent of the REFLEX sources are expected to be 
significantly contaminated by unidentified AGN. 

Figure [l] shows the spatial distribution of the REFLEX 
clusters for redshifts z < 0.35. Galactic extinction par- 
tially obscures the regions 6 h — 10 h and 16 h — 20 h . The 
cone diagrams - although averaged over a large Declina- 
tion range - illustrate the comparatively high sampling 
rates obtained with the REFLEX survey. Inhomogeneities 
in the spatial distribution of clusters on scales of the or- 
der of 100/i _1 Mpc are thus easily recognized. A detailed 
analysis of the behaviour of the mean density and of the 
topology using Minkowski functionals will be presented in 
forthcoming papers. However, the combined effect of the 
X-ray flux-limit and the steep X-ray luminosity function 
(Bohringer et al., in preparation) introduces a systematic 
dilution of the sample for larger redshifts. This is an im- 
portant difference to traditional optical cluster samples 
which are up to a certain redshift almost volume-limited 
(for given richness). In the following section the REFLEX 
data are tested for artifical number density fluctuations 
which could bias fluctuation measurements on large scales. 

3. Tests for artificial density fluctuations 

3.1. Variation of the local X-ray flux limit across the 
survey area (survey sensitivity variations) 

The local flux limit is determined by the nominal flux 
limit, the minimum number of source counts required for 
a safe detection, the local exposure time in the RASS II, 
and the local Ahi value. According to the resulting survey 
sensitivity map for > 10 source counts the nominal flux 
limit 3 x 10~ 12 erg s" 1 cm" 2 is reached on 97 percent of 
the total survey area. For > 30 source counts the fraction 
drops to 78 percent. For precise fluctuation measurements 
it is thus necessary to take into account the local survey 
sensitivity. 

In order to use as many clusters as possible for the fluc- 
tuation measurements all sources with at least 10 source 
counts in the hard band arc included. Generally, the com- 
paratively low background of the ROSAT PSPC especially 
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Fig. 3. Normalized comoving cluster number densities as a function of redshift, z, and comoving radial distance, R, computed 
with fio = 1. Vertical bars represent the formal la Poisson errors. Note the quasi-periodic density fluctuations around an 
essentially constant mean. 



in the hard band allows the detection and the character- 
ization of sources even with low source counts. In fact, 
the number of clusters with 10 to 29 source counts as ob- 
served (JVcj = 26) and as predicted from the subsample 
of the clusters with at least 30 source counts (N c \ = 37) 
suggests a formal incompleteness of 11 ± 5 clusters (ler 
Poisson error, no cosmic variance) for the subsample with 
at least 10 source counts. Assuming that the subsample 
with at least 30 source counts is complete this gives a for- 
mal overall incompleteness smaller than 3 percent. The 
corresponding local incompletenesses are expected to be 
highest in the areas where the ROSAT satellite passed 
the radiation belts in the South Atlantic Anomaly of the 
Earth's magnetic held. 

Random samples are used for the power spectrum anal- 
ysis giving Monte-Carlo estimates of the actual REFLEX 
survey windows (Sect.f|). They can also be used to test 
the quality of the survey selection model. In the follow- 
ing we describe their construction. The sensitivity map 
is computed for approximately 1° x 1° tiles covering the 
complete sky area < 2.5 deg Declination. Each of the re- 
sulting 21529 local selection functions, <fi(r), gives the 
fraction of the X-ray luminosity function at the comov- 
ing distance r, and thus the number of expected clusters, 
AN(r) = n<p(r)AV(r), down to the local flux limit of 
the given tile, 5i; m (a, 5), assuming complete randomness. 
Here, n is the mean comoving cluster number density, 
AV(r) the comoving volume element at r, and for the 
given angular coordinate (a, 5) of the tile, 



$(i X ) dLy 



where L x ln is the minimum X-ray luminosity of the sam- 
ple. For the X-ray luminosity function, $(Ix), we plug 
in the empirical estimate of the global REFLEX luminos- 
ity function as determined in Bohringer et al. (in prepa- 
ration). The shape of this function can be described by 
a Schechter function with the characteristic luminosity 
£* = 6.04 x 10 44 ergs _1 (h — 0.5), and the faint-end slope 
a = —1.61 (for a minimum of 10 source counts, no de- 
convolution of measurement error, see also de Grandi ct 
al. 1999). The transformation of the cluster restframe lu- 
minosities into the observer restframe fluxes corresponds 
to the reverse of the S — > Lx transformation described 
in Sect.|^. This prescription gives a good representation 
of the redshift histogram for L x ln = 1.0 x 10 43 ergs _1 
(comparable to the luminosity of bright Hickson groups). 



/ L ~ n <S>(L x )dL x 



(1) 



Figure U compares the observed cluster surface num- 
ber densities as a function of galactic coordinates with 
the surface number densities obtained from Monte-Carlo 
simulations of a random distribution of clusters in the RE- 
FLEX survey area with at least 10 X-ray source counts, 
and modulated by the local variation of the satellite expo- 
sure time and galactic Ahi (survey sensitivity map). The 
overall agreement is encouraging. The good statistical co- 
incidence between observed and expected cluster counts 
close to the &n = ±20° survey boundaries suggests that 
the effects of galactic extinction are well represented in the 
survey selection model. The remaining local deviations are 
caused by large-scale clustering. 
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3. 2. Variation of the average comoving cluster number 
density along the radial direction 

To test the variation of the average cluster number density 
along the radial direction, mean densities are computed for 
different volume-limited subsamples taking into account 
the local survey sensitivity map by weighting each cluster 
with the X-ray flux S using the inverse of the fraction of 
the survey area with a flux limit below S (effective area). 
For each subsample the comoving number densities are 
normalized to their respective mean density. 

Figure [| shows the normalized comoving number den- 
sity computed along the redshift direction for comoving ra- 
dial distances of R < 400 h~ x Mpc correspond to z < 0.15. 
Maximum fluctuations of the order of 3 are found on small 
scales. They are successively smoothed out with increas- 
ing R. The quasi-periodic density variations have a wave- 
length of about 150 h^ 1 Mpc. No related feature is seen 
in the power spectrum at this scale (see Sect.||). The es- 
sentially constant mean comoving cluster density implies 
the absence of selection effects discriminating against the 
more distant clusters. Note that the REFLEX survey cov- 
ers the southern hemisphere so that a volume with a ra- 
dius of R = 400 h" 1 Mpc gives a maximum comoving scale 
length of about A = 2R = 800 h" 1 Mpc. Comoving num- 
ber densities on Giga parsec scales will be discussed in 
detail in Bohringer et al. (in preparation). 

The huge nearby underdensity centered at z = 0.03 is 
also present in the ESO Slice Project data (Vettolani et 
al. 1997) as shown in Zucca et al. (1997) and might be 
the origin of the observed deficit of 'bright' galaxies in the 
magnitude number counts as discussed in Guzzo (1997). 
The large overdensity region at z = 0.05 is partially caused 
by the Shapley concentration (Fig.|l|- right cone, see also 
Scaramella et al. 1989, and Bardelli et al. 1997) and by 
some isolated nearby structures located at that distance 
in the direction of the South Galactic Pole (Fig.[l] - left 
cone). 

3.3. Flux- dependent incompletenesses 

Flux-dependent incompletenesses might also lead to sys- 
tematic errors in the fluctuation measurements. This sec- 
tion investigates the presence of this type of incomplete- 
ness and its relation to cosmic variance. 

For the REFLEX flux range the shape of the cumu- 
lative cluster number counts as a function of X-ray flux 
is mainly sensitive to flux-dependent incompleteness and 
to the K-correction, weakly dependent on evolutionary 
effects, and almost independent of the shape of a non- 
evolving X-ray luminosity function (completely indepen- 
dent for an Euclidean space), the chosen cosmological 
background model, and the type of dark matter used in the 
simulations. The comparison of the slopes of observed and 
simulated distributions provides a robust though model- 
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Fig. 4. Cumulative distributions as a function of X-ray flux 
for a REFLEX subsample (thick continuous line) and for 10 
simulated OCDM samples (thin continuous, broken, dotted, 
dashed lines) convolved with the REFLEX survey sensitivity 
and normalized to the same number of clusters. The large scat- 
ter of the number counts at high fluxes is significantly above 
the formal Poisson expectation and reflects the effects of cos- 
mic variance. 

dependent measure of the relative incompleteness of a sur- 
vey (the N-body simulations are described in Sect. |5.2| ). 

The individual cumulative flux-number counts ob- 
tained with 10 statistically independent simulations are 
shown in Fig.|]. Cosmic variance modulates the simu- 
lated cluster counts especially for X-ray fluxes S > 5.0 x 
10 -12 erg s _1 cm" 2 yielding slopes between —1.2 and —1.6. 
The fluctuations are caused by the large-scale variations 
of comoving cluster number density at small redshifts sim- 
ilar to those shown in Figs.[l] and |[ At fainter fluxes the 
fluctuations decrease and the slopes of the cumulative dis- 
tributions converge to values of about —1.3 (note that the 
plotted cumulative distributions still contain the effects of 
the effective survey area) which is close to the observed 
slope of —1.35 (Bohringer et al. 2000a). At this limit the 
REFLEX sample appears to be deep and large enough so 
that the resulting number counts should be regarded as 
statistically representative for the local Universe and not 
dominated by chance fluctuations. The similarity of ob- 
served and simulated slopes suggests a high overall com- 
pleteness of REFLEX. 

As a second measure of the overall sample incomplete- 
ness the V/Vmax test (e.g., Schmidt 1968, Avni & Bahcall 
1980) is applied as a function of the flux limit. Fig. ||shows 
the averaged V/V m &x values for different X-ray flux limits. 
Towards fainter flux limits the scatter decreases because 
sample sizes and volumes increase. At the nominal flux 
limit, the mean V/V max value is 0.512 ± 0.014 where the 
formal error does not include fluctuations caused by large- 
scale clustering. We take this convergence to the ideal case 
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Fig. 5. Mean V/V max values for REFLEX clusters as a func- 
tion of X-ray flux. Error bars are the formal ler Poisson errors 
(no cosmic variance). To illustrate the sensitivity of the test 
we have computed V/V max also below the nominal flux limit of 
3. Ox 10 -12 erg s _1 cm -2 where no REFLEX clusters are present 
to mimic a simple kind of incompleteness. Note that the test 
is performed with all clusters brighter than a given flux limit 



(0.1-2.4 keV) 



introducing a statistical dependency of the aver- 



aged VyVmax values for different S. 



< V/V max >= 0.5 for a non-expanding Euclidian universe 
as a clear sign that at the nominal flux limit the REFLEX 
survey volume and sample size is large enough to cover a 
representative part of the local Universe with a high sam- 
ple completeness and a small sample variance. 

To summarize, although it is not the basic aim of the 
present investigation to assess the absolute completeness 
and statistical representativeness of the REFLEX sample 
(see Bohringer et al. 2000a), several indications are given 
that the REFLEX survey is large enough so that in gen- 
eral the values of statistical quantities derived from the 
sample are expected to be not dominated by the effect 
of the limited REFLEX survey volume (e.g., Figs.||,||), 
and should thus give a useful characterization of the local 
Universe. The fluctuation measurements investigated here 
will not be dominated by survey incompleteness (Fig.[|) or 
other artifical large-scale variations out to radial distances 
of 400h- 1 Mpc (Fig.|). 

For the following power spectrum analyses we use dif- 
ferent subsamples which are either flux-limited (abbrevi- 
ated by F) or volume-limited (abbreviated by L). Note 
that the flux limit of the F subsamples is the nominal 
flux limit of REFLEX and that most of the F subsam- 
ples are also restricted to different volumes smaller than 
the total survey volume (see below). The characteristics 
of the subsamples are given in Tab. 1. The F0 sample con- 
tains all clusters with Lx > 10 43 ergs _1 , (h = 0.5), and 
source counts > 10. It serves as a reference sample from 
which the following subsamples are derived. The subsam- 
ples F300 to F800 differ by the chosen box length, L, used 
for the computation of the Fourier transforms, varying 
between L = SOO/i^Mpc and L = 800 hr 1 Mpc. With 
these subsamples volume-dependent effects are tested. 



The volume-limited subsamples L050 and L120 have lumi- 
nosity L x > 0.5 x lO^ergs" 1 and L x > 1.2 x lO^ergs" 1 
(h = 0.5), respectively, and are used to analyse the ampli- 
tude and shape of P{k) for clusters with different masses. 
For the given flux limit, subsamples with a lower X-ray 
luminosity cut as used in L050 are surely fluctuation- 
dominated and can thus not be regarded as statistically 
representative. L120 has the largest sample size attainable 
for volume-limited REFLEX subsamples. Tab. 1 gives co- 
moving cluster number densities and mean cluster-cluster 
distances only for the volume-limited subsamples because 
of the strong dilution of the flux-limited subsamples and 
the corresponding large change of these quantities with 
increasing redshift. 

4. Spectral analyses 

4-.1. Formal background 

In the following the spatial distribution of clusters is re- 
garded as a realisation of a formal point process. The 
corresponding Fourier transforms are well-defined in the 
strict mathematical sense if the related count measures 
are approximated by suitably smoothed versions, allowing 
the application of the classic Bochner-Khinchin theorem 
(e.g., Shiryaev 1995, p. 287) also for point processes. The 
subsequent definition of the classical Bartlett or power 
spectrum of point processes via the Fourier transform of 
reduced second-order stationary random measures (Ripley 
1977), which are closely related to the two-point (spatial) 
correlation function, does not cause any greater difficul- 
ties. More details can be found in, e.g., Daley & Vere- Jones 
(1988, Chap. 11). 

4-2. Spectral estimators 

Problems arise to find unbiased spectral estimators with 
small variance and no correlations between power spec- 
tral densities obtained at different wavenumbers k. As an 
example, naive estimators of the general form (statistical 
estimates are indicated by the hat symbol) 



P(k) 



1 



N 



N ^ 

i=i 



(2) 



where N is the number of points which are located at 
the comoving positions T*j, and Wk the discrete Fourier 
transform of the survey window, are basically applied in all 
investigations mentioned in Sect. |l[ Even for cubic survey 
volumes it leads after the subtraction of the shot noise to 
the expectations (abbreviated by the letter E) 



E{P(k)} = / F(fc-fc')P(fc')d 3 fc' . 



(3) 



where F(-) is known (in the one-dimensional case) as 
the Fejer's or Dirichlet's kernel (Percival & Walten 1993, 
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Tab. 1. REFLEX flux-limited (F) and volume-limited (L) subsamples used for the power spectral analyses of the clusters with 
> 10 source counts and the nominal flux limit 3.0 x 10~ 12 ergs -1 cm~ 2 , and Lx > 10 43 ergs _1 . The X-ray luminositities Lx 
are given for h — 0.5, box length L used for the Fourier transformation, averaged comoving cluster number densities n, and 
mean cluster-cluster distances, s, in units of h. Ncl is the number of clusters in the subsample and z the redshift. Fluxes and 
luminosities are given for the ROSAT energy band (0.1 — 2.4) keV. 



Sample 


Lx > 
[10 44 ergs -1 ] 


z < 




L 

[h,- 1 Mpc] 


n 

[h 3 Mpc-' 3 ] 


s 

[h,- 1 Mpc] 


F0 


0.1 


0.460 


428 








F300 


0.1 


0.460 


133 


300 






F400 


0.1 


0.460 


188 


400 






F500 


0.1 


0.460 


248 


500 






F600 


0.1 


0.460 


292 


600 






F700 


0.1 


0.460 


326 


700 






F800 


0.1 


0.460 


341 


800 






L050 


0.5 


0.063 


75 


400 


9.0403 x 10~ 6 


48.0 


L120 


1.2 


0.093 


96 


400 


3.8312 x 10~ 6 


63.9 



Chap. 6). The resulting systematic distortions of P(k) 
caused by the sidelobes of F increase with the dynamic 
range of P(k) and for small data volumes. Tapering is 
one method to reduce this type of leakage (Blackman & 
Tukey 1958, p. 93) but would increase the variance of P(k) 
as well. Another method is to estimate the power spec- 
tral densities only for those k values where the Fourier 
transforms of F(-) are almost zero, namely at the mul- 
tiples of the fundamental mode, fco = 2n/L, where L is 
the length of the Fourier box. In order to increase the 
signal-to-noise the power spectral densities are averaged 
over shells in fc-space with the thickness Ak — fc , cen- 
tered on the multiples of the fundamental mode. In a sim- 
ilar way smoothing of P(k) with, e.g., Bartlett, Parzen, or 
other standard spectral smoothing windows as described 
in reference books on Fast Fourier transform would re- 
duce the variance of P(k). However, as shown in Percival 
& Walten (1993), either type of smoothing is critical be- 
cause especially the central lobe of the smoothing window 
introduces a bias ~ 7 2 9 which is proportional to the 
squared bandwidth, 7 2 (a reasonable estimate of 7 is given 
by the fundamental mode fco), and to the local curvature 
of P(k). 

The leakage introduced by the survey window increases 
even further for asymmetric survey volumes because in 
this case a unique fundamental mode does not exist. For 
almost symmetric windows the effects are small and might 
be corrected using the formulae given in Peacock & Nichol- 
son (1991) and Lin et al. (1996). For highly asymmetric 
windows the whole concept of plane wave approximation 
fails. In this case the deconvolution of the survey window 
function becomes unreliable below a certain wavenumber, 
and the best solution is to resort to survey- and clustering- 
specific eigenfunctions as those provided by the Karhunen- 
Loeve transform (Vogeley & Szalay 1996). Moreover, the 
survey volume under consideration might not be large 
enough to cover a representative part of the Universe so 



that the resulting 'cosmic variance' adds to the technical 
effects described above. 

Here, for the determination of the power spectrum, two 
methods are compared. The first method uses the estima- 
tor (Schuecker et al. 1996a, b) 



P(k) 



V 



\S k \ 2 - D 



(4) 



1*1 



where the fluctuation amplitudes are corrected for the ef- 
fects of the survey window by 



1 



N 



1 



i=i 

M 



E, 



The estimator of the discreteness noise is 



D = 



EtiOX^))- 1 



(5) 



(6) 



The squared differences of the discrete Fourier transforms 
of the observed (inhomogeneous) and of the random distri- 
butions, both corrected for shot noise, are averaged over 
different directions and weighted by (1 — iWfcl 2 ) reduc- 
ing the effects of the errors in the mean number density 
(Peacock & Nicholson 1991). The power spectral densities 
must be normalized by the volume V used to compute the 
Fourier transforms and by the total power of the Fourier 
transformed survey window. Whereas the number of ob- 
served objects N is fixed by the sample, the number of 
points used for the random sample M should be large 
enough so that their shot noise contributions can be sub- 
tracted with high accuracy. Both the observed and the 
random samples have the same position-dependent selec- 
tion function, <p(r). 
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The second method to determine the power spectrum 
averages the fluctuation power over Nk modes per k shell 
(Feldman, Kaiser & Peacock 1994), 



P{k) 



1 



(7) 



where the window-corrected Fourier-transformed density 
contrasts are given in a similar way as before, 



M 



3=1 



(8) 



i=i 



The total shot noise is estimated by 

M 



V = a(l + a 



3=1 



-h-Tj 



(9) 



where a = N/M . For Gaussian fluctuations the weights 
w(r) = ! + Tt (r) p(fc) minimize the variance of the esti- 
mator, however, they require the a priori knowledge of 
P(k), that is, the quantity one wants to measure, in ad- 
dition to a fair estimate of the mean density, n(r). Rea- 
sonable results are attainable if n(r) is estimated by the 
observed luminosity function or by smoothed empirical 
z histograms (Sect. |3~l| ) and the sensitivity map of the 
survey, and if P(k) is approximated by a constant power 
spectrum, P(k) = Pq = const. 



5. Test of the spectral analyses 

5.1. General tests 

The first test concerns the choice of the spectral estima- 
tor used for the analyses of the REFLEX data. Fig. || 
compares the estimates obtained with eqs. (^) and (R). 
The power spectral densities are computed for a flux- 
limited REFLEX subsample in a cubic box with a length 
of L — 400/i -1 Mpc using a standard FFT algorithm 
on a 128 3 grid for N = 188 REFLEX clusters and for 
M = 2.0 x 10 6 random particles. The differences between 
the power spectral densities obtained with eqs. (Q) and (^) 
and the differences between the power spectra obtained 
with (|) for different P are small compa red to the errors 
introduced by the sample itself (see Sect. 5.2). We choose 
(Q) for the spectral analyses because the exploration of 
the REFLEX data should start with a minimum of pre- 
assumptions about P(k). Moreover, the REFLEX survey 
volume is comparatively symmetric so that in addition 
to the window correction term in eq. (Q) no specific de- 
convolutions are performed. The remaining effects of the 
window functions are checked using the results obtained 
with N-body simulations (see Sect. 5.2). 

To test the robustness of the method applied for the 
computat ion of the radial parts of the random samples 
(see Sect. 3A) the empirical z histogram is determined for 



500 



X[Mpc/h] 
100 50 



10 6 



10 5 



o 
a 
s 

" 10 4 



1000 r 



100 



-' 1 1 

— 


1 1 


1 


, , , | 


1 1 ' 

— 


: 


® 




















X 






xx*S$p» O" 

X 

X X * " 


r • STD 








x • 


: o fkp 


P 10x10" 


(Mpc/h) 3 






Z a FKP 


P„= 3x10* 


(Mpc/h) 3 






- x FKP 
■ 1 


P = 1x10* 


(Mpc/h) 3 
, 1 , 


, , 1 





0.01 



0.05 0.1 
k [h/Mpc] 



0.5 



Fig. 6. Power spectral densities obtained with eq.^j ('stan- 
dard method' STD) and with eq.^ ('Feldman, Kaiser, Peacock 
method' FKP) for a flux-limited REFLEX subsample with 
N = 188 clusters within a cubic volume V = (400 h' 1 Mpc) 3 . 
The radial parts of the random samples used to estimate the 
survey window are computed with smoothed empirical z his- 
tograms. 



different flux limits and smoothed with the biweight ker- 
nel (corrected for edge effects) using the standard devia- 
tion o z = 0.03 to reduce the large-scale fluctuations. The 
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Fig. 7. Fluctuation power spectral densities, P(k), as a func- 
tion of comoving wavenumber, k, corresponding to the wave- 
length A = 2Tv/k. Compared are REFLEX P{k) obtained with 
two different methods to compute the radial part of the refer- 
ence random sample. Filled symbols are obtained with random 
samples based on the X-ray luminosity function (XLF), open 
symbols are based on smoothed redshift histograms of observed 
samples (SMOOTH). 
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Fig. 8. REFLEX (continuous) and simulated (dashed-dotted) 
cluster redshift histograms. 

filtered redshift distributions give an alternative represen- 
tation of the radial selection functions (after proper nor- 
malization with the comoving volume elements and the 
survey sensitivity map). The local redshift distribution of 
the random sample as well as the local radial selection 
function is then estimated by the Monte-Carlo method. 

As an example, the power spectral densities shown as 
open symbols in Fig.^ are computed with random sam- 
ples based on smoothed empirical z distributions, the filled 
symbols with random samples based on the REFLEX X- 
ray luminosity function. It is seen that the power spectral 
densities obtained with the smoothing method are system- 
atically smaller up to factors reaching 1.6 at the largest 
scales. The differences at small scales are mainly caused by 
the poor sampling of density waves by the REFLEX clus- 
ters. We regard the luminosity function method to be more 
reliable, especially on large scales (and for small sample 
sizes): smoothing out all fluctuations is almost impossi- 
ble, especially on large scales, so that the resulting spec- 
tra have systematically smaller amplitudes as illustrated 
by Fig. 0. In the following all REFLEX power spectra ex- 
cluding those shown in Fig.|^ are obtained by using the 
luminosity function to compute the radial part of the ran- 
dom samples. 

5.2. N-body simulations 

Systematic and random errors of P(k) are investigated 
using a set of statistically independent cluster distribu- 
tions obtained from realistic N-body simulations, trans- 
formed into redshift space, and modified according to the 
REFLEX survey selection as summarized by the survey 
sensitivity map. In the following a brief overview of some 
technical aspects of the simulations are given. A more de- 
tailed description will be presented in the second paper on 
the REFLEX power spectrum. 

The simulations are performed using a standard PM 
code (Hockney & Eastwood 1988) with 256 3 particles in a 
(500 hr 1 Mpc) 3 box on a 512 3 grid giving the force resolu- 
tion Ax » l/i _1 Mpc. Ten OCDM models are simulated 



A [h-'Mpc] 
500 100 50 10 




0.01 0.05 0.1 0.5 

k [h Mpc" 1 ] 



Fig. 9. Power spectra for simulated (OCDM) flux-limited sub- 
samples in different volumes. Filled symbols give the average 
power spectral densities obtained by imposing the REFLEX 
survey conditions, continuous lines the average power spectral 
densities obtained for all-sky cluster surveys with uniform sur- 
vey sensitivities and no galactic extinction, but with the same 
X-ray luminosity function as the corresponding REFLEX sub- 
samples (error bars omitted). The error bars are the la stan- 
dard deviations of the power spectral densities obtained with 
10 REFLEX-like simulations. Error bars exceeding the plotted 
P(k) range are shown by arrows. The size of the Fourier box 
and the average number of simulated clusters located in the 
box is given in the lower left of each panel. 



with the parameters h — 0.60, cosmic density parameter of 
matter, f2o — 0.40, cosmological constant J7a = 0, cosmic 
density parameter of baryons, f2b = 0.05 (this corresponds 
to an estimate of Burbles & Tytler 1998) and as — 0.80. 
The transfer function was calculated with the Boltzmann 
code CMBFAST of Seljak & Zaldamiaga (1996). The nor- 
malization is so as to provide the correct cluster abun- 
dance satisfying both the relation given in Eke, Cole & 
Frenk (1996) and in Viana & Liddle (1996). We chose this 
model because it gives a good representation of the RE- 
FLEX data and thus realistic error estimates. However, 
any other model with a similar power spectrum could do 
the job as well. The mass resolution is 8.4 x 10 11 h~ Y M©. 
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Tab. 2: REFLEX power spectral densities obtained for different REFLEX subsamples, indicated by the index (see Tab. 1) of 
the power spectral densities, P. The errors are the formal la standard deviations as adapted from 10 OCDM simulations. 
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Each simulation starts at the redshift z = 50 (initial per- 
turbations imposed on the 'glass-like' initial load using the 
Zel'dovich approximation) and ends after 245 time steps 
(increment of the scale factor Aa = 0.004). Several repli- 
cants of the same simulation are combined using periodic 
boundary conditions to compare results on larger scales. 
However, only for scales < 500 h^ 1 Mpc statistically inde- 
pendent measurements can be obtained. 

For the identification of the clusters the friend-of-friend 
method (Davis et al. 1985) is used with the linking pa- 
rameter b = 0.16 to pick up virialized structures. The to- 
tal cluster masses are computed within the radius where 
the average density ratio is /o(r5oo)/p C riticai = 500 using 
only those clusters with at least 10 particles. The masses 
are transformed into luminosities with the empirical mass 
- X-ray luminosity relation for r$oo from (Reiprich & 
Bohringer 1999), 



h- 1 Mr. 



= 2.52 x 10 



14 



10 44 h-ievgs- 1 



0.81 



(10) 



The individual cluster masses used to derive (|10|) show 
a la scatter of about 50 percent. The r$oo radius and 
the OCDM model parameters mentioned above lead to 
realistic spatial cluster distributions, especially redshift 
histograms (Fig.^) and sample sizes, and thus to realis- 
tic error estimates for P(k). Whereas in the present in- 



vestigation the model parameters taken from the litera- 
ture are not changed, studies ar e in preparation using the 
'standard' radius r 2 oo (see Sect. 7.1 ) and different types of 
structure formation models to adjust the parameter val- 
ues in order to reconcile the models with the observations. 
We apply the above equation assuming no intrinsic scat- 
ter of the mass-luminosity relation. A simulated cluster is 
rejected if its flux is below the local flux limit given by the 
survey sensitivity map (for 10 source counts). In this way 
the simulated cluster sample follows the same sensitivity 
pattern as the observed REFLEX sample. Finally, sample 
sizes are adjusted using the additional lower luminosity 
limit Lx = 3 x 10 43 ergs _1 (h = 0.5), corresponding to 
a minimum of 46 particles per cluster. This approximate 
modeling gives realistic spatial cluster distributions and is 
enough for the error estimation of P(k). 

As a brief overview, Fig.|| illustrate the similarity of 
the observed and simulated cluster samples (see also the 
cumulative flux- number counts in Fig.Q). The model pa- 
rameters are not yet fully optimized to fit the observed 
data in detail. A more quantitative comparison is given in 
Sect.0. 

Figure^] shows the power spectra obtained from (a) 
simulated data under realistic REFLEX survey conditions 
(filled symbols), and (b) simulated all-sky cluster surveys 
with uniform survey sensitivities and no obscuration due 
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Fig. 10. REFLEX power spectra of the flux-limited subsam- 
ples F300 to F500 (filled symbols, monitored by the N-body 
simulations) and F600 to F800 (open symbols) in volumes with 
box lengths between 300 and 800 /i" 1 Mpc (no corrections for 
differences in effective biasing). For reference, the spectral fits 
obtained with the phenomenological model (continuous line) 
and with the CDM-like model (dashed line) using the subsam- 
ples F300 to F800 are superposed. 



to galactic extinction (continuous lines). The error bars of 
the latter measurements are omitted. The X-ray luminos- 
ity functions of the two sets of simulations are identical so 
that it is straightforward to test whether the power spec- 
tral estimator gives the correct power spectrum: after the 
correct elimination of the effects of the REFLEX survey 
window (see eq. ||) the resulting power spectra (shape and 
amplitude) of realistic and all-sky simulations should be 
the same. The simulations correspond to the F300 to F500 
REFLEX subsamples (Tab. 1). The errors shown in Fig. || 
represent the la standard deviations obtained from a set 
of 10 different OCDM realizations. For these simulations a 
maximum of P(k) is expected at k w 0.02 h Mpc" 1 . Note 
that the shape and amplitude of the ideal power spec- 
trum can be recovered under REFLEX conditions in all 
volumes analyzed. Some extra power is seen at the funda- 
mental mode in the 400 and 500 ft." 1 Mpc results, however 
within the la range. Note that the (500 ft" 1 Mpc) 3 simula- 
tions do not give a good representation of the fundamental 
mode at k = 27r/500 ft Mpc -1 - only 3 modes are realized 
per simulation - and do not include any fluctuations on 
larger scales, so that one should take the error bars ob- 
tained at the simulation limit with caution. Nevertheless, 
the overall agreement of the power spectra obtained un- 
der REFLEX and ideal survey conditions suggests that no 
significant systematic errors of P{k) are expected. 



6. Observed power spectra 



6.1. Exploring the general shape of P(k) 

Many variants of cosmic structure formation models dis- 
cussed today predict an almost linear slope of the power 
spectrum on scales < 40 ft" 1 Mpc and a turnover into the 
primordial regime between 100 and 300 ft" 1 Mpc. To sum- 
marize our measurements in this interesting scale range, 
Fig. [To] shows the power spectral densities obtained with 
the flux-limited REFLEX subsamples F300 to F800. The 
volumes differ by a factor 19, enabling tests of possible 
volume-dependent effects (Sect.[0|). The superposed con- 
tinuous and dashed lines in this and the following figures 
of this section are always the same. Their computation 

the 



and interpretation is described in Sect. 6.2 



In 



fol- 



lowing they may serve as a mere reference to compare 
the power spectra obtained with the different REFLEX 
subsamples listed in Tab. 1. Fig.|ll| gives a more detailed 
view of the spectra obtained with the subsamples F300 
to F500 in volumes which are monitored by our N-body 
simulations. Fig. [l^ compares the spectra obtained for the 
volume-limited subsamples L050 and L120 with the spec- 
trum obtained for the flux-limited subsample F400, all 
spectra are estimated within the same Fourier volume. 
Finally, Fig.|l^ shows the combined power spectrum ob- 
tained with the subsamples F300 to F500 which we regard 
as the basic result of the REFLEX power spectrum anal- 
yses. The values of the power spectral densities obtained 
with the subsamples F300-F500, L050, and L120 with the 
errors estimated with the N-body simulations are given in 
Tab. 2. In the following a few more detailed remarks are 
given. 

Figure [l0| shows the superposition of the power spectra 
obtained with the flux-limited subsamples F300 to F800 
in the comoving volumes ranging from (300 ft" 1 Mpc) 3 to 
(800 ft" 1 Mpc) 3 . The data are not corrected for sample-to- 
sample variations of the effective biasing (see Sect. [7]) so 
that the effective variance among the estimates is possi- 
bly smaller than that displayed by the figure. The point 
distribution outlines a corridor which can be separated 
into three parts. For k > 0.1 h Mpc" 1 the power spec- 
tral densities decrease approximately as fc" 2 . Between 
0.02 < k < 0.1 h Mpc -1 the spectra bend into a flat 
distribution. The N-body simulations give la standard 
deviations between 30 and 80 percent (including cosmic 
variance) in this scale range as shown in Figs.[ll], [l2[ 
||. For k < 0.02 h Mpc" 1 a second maximum is seen at 
k w 0.01 h Mpc" 1 . We did not perform N-body simula- 
tions for such large scales. However, the delete-d jackknifc 
resampling method (a variant of the boostrap method 
where the creation of artifical point pairs is avoided; see, 
e.g., Efron & Tibshirani 1993, see also the critical remarks 
on the use of the bootstrap method in point process statis- 
tics given in Snethlage 2000) gives la error estimates of 
the order of 80 percent (cosmic variance not included). 
The detection of the second maximum in the power spec- 
trum on such large scales, if real, would have very im- 
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Fig. 11. REFLEX power spectra of the flux- limited subsam- 
ples F300 to F500. The box lengths and the number of clusters 
used for the power spectrum estimation are given in each panel 
in the lower left. The bars represent the la errors adapted from 
N-body simulations. The fits of a linear power spectrum model 
(dashed lines) and of the phenomenological model (continuous 
line) using the subsamples F300 to F800 are superposed. 



portant implications on current structure formation mod- 
els. However, as pointed out in the Introduction, measure- 
ments on such large scales are easily biased by very small 
systematic errors of the survey detection model. We post- 
pone a detailed study of this very questionable feature to a 
subsequent paper. The present investigation concentrates 
more conservatively on the range 0.013 < k < 0.4 h Mpc -1 
which is found to be free from any significant artifical fluc- 
tuations (Sect.||), which can be easily monitored by the 
available N-body simulations, and which contains density 
waves well sampled by the REFLEX clusters. 

Individual spectra obtained with the three flux-limited 
subsamples F300 to F500 are shown in Fig. O, now includ- 
ing the ltr errors adapted from the N-body simulations. 
Whereas the spectra obtained with F300 and F400 (upper 
and middle panel) show a maximum at k « 0.03 h Mpc -1 , 
the F500 data (lower panel) suggest only a flattening of 
the spectral densities. Especially the power spectral den- 
sity obtained at the fundamental mode seems to indicate 
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Fig. 12. REFLEX power spectra of the volume-limited sub- 
samples L120 (upper panel), and L050 (middle panel), and of 
the flux- limited subsample F400 (lower panel). LI 20 contains 
clusters with a brighter lower X-ray luminosity compared to 
L050 (and F400). The bars represent the la errors adapted 
from N-body simulations. For reference the fits obtained with 
the phenomenological model (continuous lines) and with the 
CDM-like model (dashed lines) using the subsamples F300 to 
F800 are superposed. The amplitudes of P(k) increase with in- 
creasing lower luminosity limit as expected by standard biasing 
schemes. 



a still rising power spectrum for smaller k values. A simi- 
lar effect is seen in the simulations (see Fig.^) suggesting 
a statistically not very significant but noticable leakage of 
fluctuation power especially from the second to the first 
fundamental mode. The reference to Fig. |l0| reveals that 
the fundamental mode of F500 is already part of the sec- 
ond probably not real maximum in the power spectrum at 
k w 0.01 h Mpc -1 . Hence the fundamental mode of F500 
should not necessarily get the highest weight in the eval- 
uation of the maximum of the power spectrum on smaller 
scales. We test the possibility that the location of the max- 
imum increases with volume but could not find any sys- 



tematic effect (see Sect. 6.2). 

The spectra shown in Fig. ^2] obtained with the 
volume-limited subsamples L050 and L120 (upper and 
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Fig. 13. Combined REFLEX power spectrum obtained with the flux-limited subsamples F300 (open squares), F400 (filled 
hexagons), F500 (open hexagons), and their standard la deviations adapted from N-body simulations. Not shown are the 
power spectral densities with wavenumbers k > 0.4 ft Mpc -1 because the corresponding density waves are only sparsely sampled 
by REFLEX. The continuous line is a fit of the phenomenological model, the dashed line the CDM-like model fit using the 
flux-limited subsamples F300 to F800. The power spectral density at 500 ft -1 Mpc is part of the extra power detected at 
k « 0.01 h Mpc -1 and might already be biased. 



middle panel) show a broad maximum at k » 
0.03 /iMpc -1 . A weak indication is found for a positive 
slope on larger scales. The second maximum of the power 
spectrum seen in Fig.[l(] is not sampled by L050 and L120 
because their sample volumes do not reach such large 
scales. The Fourier volumes are therefore restricted in 
both cases to (400 h Mpc -1 ) 3 . For comparison the lower 
panel shows the power spectrum obtained with the flux- 
limited subsample F400 estimated within the same volume 
as used for L050 and L120. In general, the overall shapes 
of the spectra obtained with the volume- and flux-limited 
subsamples are found to be similar, although minor dif- 
ferences might be seen on smaller scales (see below). The 
three spectra also show that the amplitude increases with 
increasing lower X-ray luminosity of the subsample. How- 
ever, larger sample sizes are needed to confirm the effect. 

To summarize, basically all REFLEX spectra are con- 
sistent with a broad maximum of the cluster power spec- 
trum at comoving wavenumbers around k w 0.03 h Mpc -1 
corresponding to wavelengths of about 200 h" 1 Mpc. A 
second maximum is found at k = 0.01 h~ l Mpc corre- 
sponding to 600 ft- -1 Mpc, but appears questionable (see 
Sect.H). These findings are summarized in Fig. 13, show- 



ing the combined spectra obtained with the subsamples 
F300 to F500, and illustrating the stability of the results 
obtained within different volumes. We regard this as a 
representative REFLEX power spectrum. 



6.2. First cosmological implications 

In the following we want to characterize the overall shape 
of the observed power spectra as well as specific spec- 
tral features like the location of the maximum and the 
local slope of P(k) in specific k ranges, restricting the dis- 
cussion mainly to the conservative k range 0.013 < k < 
0.4 h Mpc -1 mentioned above. This will enable us to de- 
rive our first cosmological implications. 

The REFLEX power spectral densities shown in the 
last sections are sampled strictly following the rules of 
standard Fourier analysis. As a consequence we have to 
work with uncomfortably large but statistically almost in- 
dependent k bins which complicates the analyses even of 
the maximum of P(k) in the conservative k range. To im- 
prove the 'eye ball' estimates of the locati on o f the max- 
imum of the power spectra given in Sect. |6. l| and to get 
a handle of the expected errors, the spectra are param- 
eterized in two different ways. The first method applies 
a purely phenomenological fitting function which gives an 
almost model-independent description of the data (see also 
Peacock 1999, p. 530): 



P(k) 



A k~' 11 k il ~ 3 



(11) 
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The location of the maximum of the power spectrum is 
i 

-n\ " _Tls 



h — h 



(12) 



The slope on large scales, n — i\ + i 2 — 3, is set to '1' 
because no statistically reliable information is attainable 
from the REFLEX data in this scale range. The slope on 
small scales is n s = ii — 3. The characteristic scale, fc ma x, 
is comparable to the wavenumber corresponding to the 
horizon length at the epoch of matter-radiation equality, 
fc cq = 0.195 Oo^Mpc" 1 (see Peebles 1993, p. 164, and 
Peacock 1999, p. 459), and thus yields an estimate of the 
cosmic mass density (assuming that 3 relativistic neutrino 
families are left over from high redshift, and that neutrino 
masses are small compared to the temperature of the cos- 
mic microwave background radiation), 



J7n h 



3.13- 



fen 



hMpc 



(13) 



In contrast to the first method which allows a variable 
slope at small scales and a narrow maximum, the second 
method is less flexible, but physically better defined. The 
fitting function is based on the CDM linear transfer func- 
tion, T(fc), as given in Bardeen et al. (1986), where the 
power spectrum, P{k) = AkT 2 (k), again is assumed to 
have n = 1 on very large scales. The shape of the power 
spectrum is characterized by the shape parameter, T, de- 
fined in the standard way (eq.|l^), giving the approximate 
relation between T and fc max 



0.114r/iMpc~ 



(14) 



For CDM-like models with low baryon density, T is mainly 
determined by f2o, h, and f2b (see eq.|l9|). This equation 
offers another way to approximate the cosmic mass density 



flnh w r. 



(15) 



A standard SIMPLEX \ 2 minimization method is applied 
separately to the spectra obtained with the subsamples 
F300 to F800 to perform numerical fits from which the 
values of fc max and T are deduced. This assumes that the 
power spectral densities of each individual spectrum are 
statistically independent. For the given REFLEX survey 
window (|W4| 2 < 0.082 for all k and volumes studied), 
and for the given spacing of the k values of the measured 
P(k) data at the multiples of the fundamental mode this 
is approximately the case. The values of £; max and T ob- 
tained from the fits are independent of the volumes used to 
perform the Fourier analyses as shown in Fig. |l4|, strongly 
supporting the detection of a real maximum of P(k) in 
the given k range. Averages and their formal la standard 
deviations of fc max an d T using the subsamples F300 to 
F800 give for the two fit functions, respectively, 



fcnn.x = 0.030 ± 0.005. 



T = 0.195 ± 0.055. 



(16) 
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Fig. 14. Upper panel: values of the shape parameter, T, as a 
function of the box length, L, of the Fourier volume, obtained 
from fits of the linear CDM model to the flux limited subsam- 
ples F300 to F800 (from left to right). Lower panel: values of 
the wavenumber of the maximum of the power spectrum, fc max , 
as a function of L obtained from fits of the phenomenological 
model. The sample averages are indicated by dashed lines. 

The T estimate corresponds to fc max = 0.022 ± 0.006. 
Similar numbers are obtained when only the subsamples 
F300 to F500 are used. Note that the subsamples F300 
to F800 are statistically dependent so that the error esti- 
mates given in ( |l6| ) must be regarded as lower limits. The 
phenomenological and the CDM-like model based on these 
mean values are shown as continuous and dashed lines, re- 
spectively, in the Figs. [H)[ [ll], || || and || They both 
give a good description of the shape of all power spec- 
tra obtained with the flux- and with the volume-limited 
REFLEX subsamples. In Fig. [l^ we show this for the com- 
bined power spectrum obtained with the flux-limited sub- 
samples F300 to F500. The two methods give consistent 
results for the location of the maximum of the REFLEX 
power spectra in the range 



0.022 ±0.006 < /c max < 0.030 ± 0.005 /iMpc" 1 



(17) 



Similarily, the values of the cosmic density parameter ob- 
tained with eq. ( fl3| ) and /c max from (fT6|), and with eq. (15) 
and r from (jig) give the range 



0.15 ±0.03 < flnh < 0.20 ±0.06. 



(18) 



It is interesting to note that the spectral slopes on 
small scales of the volume-limited subsamples L050 and 
L120 estimated by fitting ( pd| ) are found to be slightly 
flatter, n s = —1.6 ± 0.4, compared to n s = —1.8 ± 0.4 
obtained with the flux-limited subsamples (ler error esti- 
mates from N-body simulations). However, the differences 
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Fig. 15. Combined REFLEX power spectrum obtained with 
the subsamples F300 to F500 (filled symbols) compared to 
the power spectrum obtained from Abell/ACO clusters (open 
hexagons) by RetzlafT et al. (1998) and from APM clusters 
(open squares) by Tadros et al. (1998). 



are statistically perhaps not very significant and might be 
attributed to cosmic variance. 

Figure |l5| compares the REFLEX power spectrum with 
the Abell/ACO (Retzlaff et al. 1998, see also Einasto 
et al. 1997) and the APM (Tadros et al. 1998) spectra. 
The respective amplitudes of the power spectra of the 
Abell/ACO and APM samples are 1.7 and 2.2 below RE- 
FLEX. This might be attributed to the different cluster lu- 
minosities contained in the samples. For k > 0.08 h Mpc -1 
the spectra give consistent slopes of approximately —1.8 
although both the REFLEX and the Abell/ACO sample 
do not show the minimum at k ~ 0.1 /iMpc -1 found with 
the APM sample. Regarding the maximum of P(k) the 
Abell/ACO data suggest a comparatively narrow peak 
at fc max = 0.05 /iMpc -1 consistent with the estimate of 
Einasto et al. (1997). Contrary to this the REFLEX spec- 
trum has a broad maximum which peaks in the range 
0.022 < fc max < 0.030 /iMpc -1 . Note that the exact eval- 
uation of the statistical significance of this difference is 
difficult to assess because the REFLEX and Abell/ACO 
power spectra are sampled in different ways. The broad 
maximum of the REFLEX spectrum appears to be more 
consistent with the APM sample if the REFLEX measure- 
ment at 500/i _1 Mpc is excluded. 

Figure |l6| compares the combined REFLEX power 
spectrum obtained with the flux-limited subsamples F300 
to F500 with the spectrum obtained with a magnitude- 
limited sample of Durham/UKST galaxies (Hoyle et al. 
1999). We chose this sample because of the comparatively 
large samples size (2501 galaxies, 1 in 3 sampling rate), 
the large volume (1450 square degrees, z < 0.1), and the 
small effects of the survey window. Recall that the upper 
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Fig. 16. Combined power spectrum obtained with the RE- 
FLEX cluster subsamples F300 to F500 (filled symbols, mea- 
surements on scales < 20/i -1 Mpc omitted) compared with 
the power spectrum of Durham/UKST galaxies (open sym- 
bols) obtained by Hoyle et al. (1999). The upper lines give the 
fit to the REFLEX cluster power spectra, the lower lines the 
same fits divided by the squared 'biasing factor' b 2 = 2.6 2 . The 
continuous line is a fit of the phenomenological model (Par), 
the dashed line the CDM-like model fit (CDM) as described in 
Pig.g 



continuous line is the fit of the phenomenological model 
to the REFLEX data, the upper dashed line the fit of the 
CDM-like model; the lower lines are the same fits shifted 
by the factor 6.8. For wavelengths 20 < A < 300 h^Mpc 
the overall shapes of the cluster and galaxy power spec- 
tra are very similar. The ratio of the linear biasing factors 
for the given REFLEX cluster subsample and the galaxy 
sample as deduced from the shift factor is b = 2.6. 

7. Comparison with CDM models 

7.1. Semi-analytic model 

To make a first comparison with cosmological models and 
an attempt to differentiate between their presently dis- 
cussed variants, an outline of a semi- analytic model is 
given for biased nonlinear power spectra in rcdshift space 
for clusters of galaxies. The model gives a good overview 
of the effects of different model parameters and is used 
to narrow the parameter ranges needed for a more de- 
tailed comparison with N-body simulations. Notice that a 
significant number of N-body simulations has to be per- 
formed for each parameter set in order to derive statistical 
meaningful error estimates which is planned for the sec- 
ond paper on the REFLEX power spectrum. The model 
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spectra are computed with parameter values taken from 
the literature and are compared with the REFLEX power 
spectra. No evolution of structures is assumed within the 
redshift range covered by the REFLEX subsamples ana- 
lyzed (z < 0.15, for an exact treatment see also Magira 
et al. 2000). The linear power spectrum P(k) ~ k n T 2 {k) 
is normalized by the standard deviation of the density 
contrast, erg, obtained with the spherical top hat filter 
function with the filter radius R = 8/i -1 Mpc. The fit- 
ting formula for the linear transfer function, T(k), from 
Bardeen et al. (1986) is applied with the shape parameter 
(Sugiyama 1995) 



r = Qq h exp 




(19) 



appropriate for models with low present baryon densities, 
fib ■ The mapping between evolved and non-evolved real- 
space power spectra can be deduced from fits to the re- 
sults obtained with N-body simulations (Hamilton et al. 
1991, Peacock & Dodds 1994). We apply the prescription 
in the form presented in Mo, Jing & Bonier (1997, and 
references given therein). The approximation of the lin- 
ear growth factor at redshift zero is taken from Carroll, 
Press & Turner (1992). It was shown that this mapping is 
almost independent of the slope of the linear power spec- 
trum. Moreover, as long as the strongly non-linear cluster- 
ing regime is excluded the formalism reproduces the power 
spectra obtained from N-body simulations quite well. 

To compute the observed or effective biasing values 
Moscardini et al. (2000, see also Matarrese et al. 1997 and 
Borgani et al. 1999) assumed a linear biasing between mat- 
ter and object number density fluctuations, a reasonable 
assumption in the linear regime. They derived an exact 
relation between the observed and the matter two-point 
spatial correlation function (their equation 7) which we 
reproduce in /c-space, ignoring any redshift-dependence of 
the correlation function. The real-space (evolved) power 
spectrum thus reads 



PreAk) = P(k) 



dz 1 N(z 1 )r~ 1 (z 1 ) 



(20) 



dz 1 dz 2 N(zi)r 1 (zi)b cS (z 1 )b eS (z 2 )r 1 (z 2 )N(z 2 ). 



where N(z)dz is the number of clusters expected in the 
redshift interval [z, z + dz], 



bett(z) 



1 



N(z) 



dM b(M) N(M, z) 



(21) 



the effective biasing at redshift z, N(M, z)dM dz the num- 
ber of clusters expected in the mass range [M, M + dM] 
and in the redshift range [z, z + dz], and r(z) the comoving 
distance of the redshift shell z. Note that N(M, z)dM = 
N(L, z)dL so that © can also be used for luminosi- 
ties once the M — L conversion is performed (see below). 



We evaluate the integrals at the redshifts and luminosi- 
ties given by the observed subsample for which the effec- 
tive biasing factor should be determined (for each clus- 
ter we compute from the luminosity the mass and the 
corresponding biasing factor and plug the result into the 
two equations given above). This circumvents the intro- 
duction of additional Press-Schechter -like models for the 
cosmic mass function which we do not intend to test in the 
present context (probably introducing some inconsistency 
between the observed and the model luminosity function), 
but guarantees a one-to-one correspondence between the 
clusters used to measure the power spectrum and the clus- 
ters used to estimate the effective biasing for the subsam- 
ple under consideration. This, however, increases the vari- 
ances but reduces the systematic errors to a minimum. 

For clusters of galaxies simple biasing schemes are ex- 
pected (Sect. In this respect the model of Mo & White 
(1996) is of special interest. They combine (a) conditional 
probability densities derived by Bond et al. (1991) for 
Gaussian random fields within the general framework of 
Markovian diffusion processes with an 'absorbing barrier' 
at the critical density contrast, with (b) gravitationally in- 
duced motions as predicted by a spherical collapse model. 
We use the fitting formula given in Sheth & Tormen (1999) 
which is found to give a better agreement with N-body on 
small scales. The critical overdensity, 5 C , is determined by 
the cosmological background model as described in Ki- 
tayama & Suto (1996). The relation between mass and 
radius is M = °°f| fl3 . 

For the M — Lx conversion the empirical relation be- 
tween the total mass M and X-ray luminosity within t 2 qq 
is used (Reiprich & Bohringer 2000): 



M 



Mr. 



4.7 x 10 1 



Ly 



10 44 h~ 2 ergs" 1 



(22) 



The systematic effects caused by using the relation ob- 
tained with r 50 o are investigated below (Fig.|l8|). We as- 
sume that the masses deduced from ( p2| ) closely resemble 
the virial masses. 

The transformation of the real-space power spectrum 
into redshift space is determined by the effects of peculiar 
velocities and redshift measurement errors. If the max- 
imum distances are large compared to (distant ob- 
server approximation) only the linear flow of the velocity 
field makes an additional contribution to the fluctuation 
field in redshift space (Kaiser 1987). On small scales the 
peculiar velocities and the redshift measurement errors of 
the clusters smooth the fluctuation field which can be de- 
scribed by a Lorentzian distribution in fc-space. The two 
effects can be integrated over the cosine, fi, of the angle 
between the normal vector of the density wave in fc-space 
and the line-of-sight, and give 



Pobs(k) = Preal(fc) / cZ/1 



(i + /V) 5 



o r (l + fcV^) 



(23) 
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(24) 



Here, op is the pairwise velocity dispersion of the dark 
matter haloes and o~ z the average error of the cluster red- 
shifts. The effects of S1a are not important (see Lahav et 
al. 1991). For the computation of o~p we concentrate on 
the effects on scales (pair separations) larger 5 ft" 1 Mpc 
and neglect correlated motions so that <7p can be approx- 
imated by (Peebles 1980, Sect. 72) 



< vr > 



(25) 



It is thus expected that the cluster motion inherit a ran- 
dom velocity from the random motion of the overall mat- 
ter distribution. The density-weighted mean square pe- 
culiar velocity is determined by the integral of the cosmic 
energy equation as derived in the BBGKY hierarchy (e.g., 
Peebles 1980, Sect. 74, Mo et al. 1997): 



< vt > 



- fl H a 



dk 



A|(fc). 



(26) 



The cluster redshift errors can be deduced from the red- 
shift errors of the individual galaxies (a g = 200 km s" 1 , 
Guzzo et al., in preparation), the number of cluster 
galaxies used to estimate the cluster redshift {N z — 5) 
and the line-of-sight cluster velocity dispersion (ctlos = 
700 kms -1 , see Zabludoff et al. 1993). Although the lat- 
ter quantity depends on the specific structure formation 
model we take the empirical estimate because in the 
present case we are only interested in a redshift error esti- 
mate where the effects of the structure formation models 
are of second-order. Following Danese, deZotti & diTullio 
(1980) the squared error of the cluster redshift is 



[1.178 a£ os + a g 2 ] 



(27) 



where the factor 1.178 results from the 68 percent confi- 
dence interval of the Student's t-distribution with 4 de- 
grees of freedom (for N z = 4, 3, 2 we computed the factors 
1.242, 1.391 and 2.057, respectively). Typical redshift er- 
rors of the REFLEX clusters of galaxies are thus expected 
to be <j z — 350 kms" 1 . 

7.2. Test of the model 



We test (g5j) and (|2J) with the available OCDM N-body 
simulations (parameters are giv en in Ta b. 3), yielding the 
mean cluster peculiar velocity \J < v \ > = 516 kins" 1 (the 
standard deviation of this quantity for different simula- 
tions is 3 percent), and for pair separations > 10 ft" 1 Mpc 
the approximately constant (within about 3 pecent) pair- 
wise cluster velocity dispersion op = 430 km s" 1 . Within 
the given errors the relation between these quantities is 
reproduced by (|25|). The simulations give a maximum 
at 3 ft." 1 Mpc of op = 500 km s" 1 and values of about 
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Fig. 17. Squared biasing factors computed with eq. (^8|) for 
the F0 subsample as a function of pair separation r. The la 
error bars include the effects of the errors introduced by the 
50 percent uncertainty of the X-ray mass-luminosity relation. 
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Fig. 18. Test of the semi-analytic model with N-body simula- 
tions. Shown is the power spectrum averaged over 10 OCDM 
N-body simulations (filled symbols) of ideal clusters samples 
(Sect.H) and their la standard deviations (error bars). The 
lines represent the power spectra obtained with the semi- 
analytic model for the same model parameters as the N-body 
simulations (fio = 0.40, Qa = 0, D. h = 0.05, h = 0.60, F = 0.20, 
n — 1, as = 0.80, i.e., cluster-normalized). Dashed-dotted line: 
linear matter power spectrum. Long-dashed line: evolved mat- 
ter power spectrum. Short-dashed lines: evolved power spec- 
trum including effective biasing. Continuous lines: last spec- 
trum transformed into redshift space. The two types of spectra 
are shown for the mass-luminosity relation obtained with r2oo 
and r 50 o- 

200 km s" 1 on smaller scales. The semi-analytic model 
neglects the small-scale dependency because the RE- 
FLEX power spectrum does not sample the correspond- 
ing k range. On the other hand using (^5|) and (^6[) the 
semi-analytic model predicts ap = 360 km s" 1 , which is 
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about 15 percent too small compared to the simulations. 
We found this approximation good enough for the real- 
rcdshift space transformation. 

An important assumption implicitely used for the 
derivation of fl20| ) is that the averaged biasing factor is in- 
dependent of pair separation, r. For flux-limited samples 
one might expect that at large r the fraction of pairs con- 
sisting preferentially of at least 1 luminous cluster could 
artifically increase the effective biasing factor. This would 
increase the measured power spectral densities at small 
k and thus steepen the slope compared to the volume- 
limited case. To test this, the number of pairs with sepa- 
ration r are weighted with the individual biasing factors 
of the pair members, yielding the average squared biasing 
factors, 



Tab. 3. Model parameters of CDM variants used for the semi- 
analytic model. 



< b 2 (r) > = 



E 



\r< \ri — Tj | <r+Ar} 



b(Mi) b{M 3 



N P (r) 



(28) 



where iVp(r) is the number of cluster pairs with separa- 
tions within [r, r + Ar]. The mass variances, <r 2 (M), used 
for the computation of 6(Mj) are derived assuming a scale- 
invariant power spectrum with the spectral index —2. 
Fig. [l?] shows the average squared biasing factors as a func- 
tion of pair separation for the REFLEX subsample FO. For 
pair separations r < 150 h~ x Mpc no scale-dependent cor- 
relations between the individual biasing factors are seen. 
For larger separations the systematic increase of the aver- 
age squared biasing factor suggests that the treatment of 
effective biasing as described above must be modified. The 
maximum pair separation, r, corresponds to the minimum 
wavenumber k k, n/2r — 0.010 /iMpc^ 1 . For wavenum- 
bers larger than this limit no systematic errors are ex- 
pected. It will be seen that the observed REFLEX power 
spectra which are compared with the biasing model do not 
reach this limit. Note that this refers only to REFLEX and 
must be re-evaluated for other surveys. 

The semi-analytic model is tested against the 10 N- 
body simulations (OCDM) of ideal cluster samples de- 
scribed in Sect.|| In Fig.[l8] the lines give the theoretical 
spectra obtained under the different model assumptions, 
the filled symbols the average power spectral densities ob- 
tained from the N-body simulations, and the error bars 
their la standard deviations. The overall agreement be- 
tween model and simulation is good enough to separate 
between different scenarios of cosmic structure formation. 
The largest ambiguity is introduced by the specific choice 
of the mass-luminosity relation. In the following the the- 
oretical spectra obtained with 7-200 are shown because the 
corresponding cluster masses are expected to give better 
estimates of the virial masses. 

7.3. Results 

As an example, in Fig. the REFLEX power spectrum 
obtained with the F400 subsample is compared with differ- 
ent variants of CDM models (the data obtained with F300 



Model Q Ma 



(J* 



SCDM 1.00 0.00 0.50 1.0 0.050 0.45 1.37 

OCDM 0.40 0.00 0.60 1.0 0.050 0.20 0.80 

ACDM 0.30 0.70 0.65 1.0 0.036 0.21 0.93 

TCDM 1.00 0.00 0.50 0.8 0.100 0.41 0.58 

tCDM 1.00 0.00 0.50 1.0 0.050 0.21 0.60 



and F500 give similar results). The values of the model 
parameters are given in Tab. 3. The standard Cold Dark 
Matter (SCDM) model with the COBE normalization as 
given in Bennett et al. (1994) is shown for reference. 
The open CDM (OCDM) model is cluster-normalized (see 
Sect.||). For the low-density flat (ACDM) model see Lid- 
die et al. (1996a,b). The tilted (TCDM) model is de- 
scribed in Moscardini et al. (2000) and the references given 
therein. The rCDM model is cluster-normalized according 
to Viana & Liddle (1996). 

The measured power spectra discriminate between the 
models, SCDM and TCDM are excluded, rCDM fits 
marginal the lower la range, the open and ACDM models 
slightly underpredict the fluctuation amplitude but within 
the la significance range. 

To test the biasing trends we changed the ACDM nor- 
malization from as = 0.93 to ag = 0.70 (similarly we 
could also change a 8 = 0.80 to a$ = 0.60 for the OCDM 
model) yielding an acceptable fit to the flux-limited RE- 
FLEX power spectrum (open symbols and continuous line 
in Fig. ^). The ACDM spectra are then computed for 
the same volume-limited subsamples as used for the de- 
termination of the empirical spectra. The increase of the 
amplitude with the increasing lower X-ray luminosity - 
although at the detection limit of REFLEX - is well re- 
produced by the model, but not the apparent flattening 
of the slope on scales < 100 /i -1 Mpc. However, the er- 
rors of the slope measurements as deduced from the sim- 
ulations are quite large so that the apparent difference 
might not be statistically significant. Moreover, neither 
the scale-independency of the effective biasing parameter 
in this range (see Fig.[l7]) nor the analyses of the OCDM 
simulations suggest such an effect. 

8. Discussion and conclusions 

The most important result of the present investigation 
is the detection of a broad maximum of the power spec- 
trum of the fluctuations of comoving number density 
of X-ray selected cluster galaxies in the range 0.022 < 
k < 0.030 /iMpc -1 (Fig.rfl). The maximum is flatter and 



peaks at a smaller wavenumber compared to optically se- 
lected cluster samples. On scales 0.02 < k < 0.4 h Mpc" 1 
the similiarity to the spectra obtained from optically se- 
lected galaxy samples is striking (Fig.[l6|). In this range 



20 



The REFLEX Power Spectrum 
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Fig. 19. Comparison of observed power spectral densities and 
expectations of variants of CDM semi-analytic models for the 
flux-limited subsample F400. The model parameters are given 
in Tab. 2. 
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Fig. 20. Comparison of observed power spectral densities and 
predictions of the ACDM semi-analytic model for the flux- 
limited subsample F400 (open hexagons) and for the volume- 
limited subsamples: stars for Lx > 0.5 x 10 44 ergs _1 (subsam- 
ple L050), filled symbols for Lx > 1.2 X 10 44 ergs -1 (subsample 
L120), h = 0.5. The A model is renormalized to erg = 0.70 to 
give a good fit to the flux-limited subsample. 



the REFLEX data rule out galaxy formation models with 
strongly nonlinear biasing schemes. 

Within the course of the exploration of the REFLEX 
data and the results of the N-body simulations we found 
that for surveys smaller than REFLEX cosmic variance 
might be more important than previous studies suggest. 



For example, the variation of the comoving cluster num- 
ber density along the redshift direction shows a huge un- 
derdense region located between z as 0.015 and 0.045 in 
the southern hemisphere where the comoving cluster num- 
ber density drops by a factor of 3 below the mean level 
(Fig.||). This complicates the determination of the local 
cluster luminosity function, at least for the less rich sys- 
tems (Bohringer et al., in preparation). Another example 
is the variation of the linear slope of the cumulative flux- 
cluster number counts between —1.6 and —1.2 as found 
in the N-body simulations (Fig.^J). We regard this as a 
warning not to draw general cosmological conclusions from 
cluster samples with a size smaller than REFLEX. 

The REFLEX data show extra fluctuation power on 
scales k « 0.01 hMpc^ 1 (Fig.^). From our simulations 
we found that artifical power spectral densities of an or- 
der of magnitude can be easily produced on 500 h~ x Mpc 
scales if, e.g., the lower X-ray luminosity limit of LyP a = 
1.0 x 10 43 ergs _1 , which is used in the present investiga- 
tion to get almost complete REFLEX subsamples, would 
be erroneously underestimated by a factor of about 1.5. 
Similarily, already on scales of 400 h~ x Mpc small changes 
in the method to estimate the radial part of the selection 
function (compare the results obtained with smoothed 
redshift distributions and X-ray luminosity functions, 
Fig. 0) change the power spectral densities by a factor 
1.6. These two examples illustrate the difficulty measur- 
ing fluctuations on scales > 400 h^ 1 Mpc which is the basic 
motivation for restricting the present investigations more 
conservatively to the small wavelength range. 

Extra fluctuation power on 800 /i~ 4 Mpc scales is also 
found for the Abell/ACO richness > 1 clusters by Miller 
& Batuski (2000). In addition to the fact that they over- 
sample the cluster power spectrum which mimic a more 
significant effect than the data can provide, it is difficult 
to understand how gradients in comoving cluster number 
density by a factor of 2, corrected with crude step- like ra- 
dial selection functions, and the neglection of any correc- 
tions for galactic extinction can lead to precise fluctuation 
measurements at 800 ft. -1 Mpc. It is surely insufficient to 
use cluster quadrant counts showing a scatter of 16 per- 
cent to justify fluctuation measurements aiming to detect 
fluctuations below the 1 percent level. 

The REFLEX power spectra do not show any indica- 
tion for a narrow peak at k = 0.05 h Mpc - . The report 
of such a feature in the power spectrum of Abell/ACO 
clusters and the interpretation as evidence for a regu- 
lar distribution of galaxy clusters with a periodicity of 
120 /i~ 4 Mpc by Einasto et al. (1997) implies substan- 
tial difficulties for current models of structure formation. 
Retzlaff et al. (1998) who have found a similar but less 
peaked feature in the Abell/ACO cluster P(k) used a 
large set of N-body simulations to demonstrate the po- 
tential importance of cosmic variance in this context. The 
discrepancy between REFLEX and Abell/ACO cluster re- 
sults might be attributed to the additional 35 percent 
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non-Abell/ACO/Supplement clusters included in the RE- 
FLEX catalogue. Unfortunately, the subtle selection ef- 
fects imposed by optical cluster selection (Sect.[l]) makes 
a quantitative discussion of this point almost impossible. 
In any case, due to current sample depths, cluster power 
spectrum analyses are restricted in general to volumes 
< (500 h^ 1 Mpc) 3 , and this imposes a spectral resolution 
Ak = 0.013 (fundamental mode) at best. Therefore, a sig- 
nificant detection of a feature such as a peak of width 
Ak « 0.02 is arguable at all. 

The REFLEX spectra are compared with semi- 
analytic models describing the biased nonlinear power 
spectrum in redshift space. Most of the equations applied 
are calibrated with N-body simulations. We found that 
structure formation models with a low cosmic mass den- 
sity (OCDM, ACDM) give the best representation of the 
REFLEX data (Fig. ^9j) . Although the models could repro- 
duce the observed changes of the amplitudes with samples 
of different luminosities, we regard the results are tenda- 
tively. Larger sample sizes are necessary to confirm this 
finding. 
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